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Abstract: 

New fundamental particles, charged under new gauge groups and only weakly coupled to 
the standard sector, could exist at fairly low energy scales. In this article we study a 
selection of such models, where the secluded group either contains a softly broken U{\) 
or an unbroken SU{N). In the Abelian case new 7„ gauge bosons can be radiated off 
and decay back into visible particles. In the non-Abelian case there will not only be a 
cascade in the hidden sector, but also hadronization into new vr^, and /9„ mesons that can 
decay back. This framework is developed to be applicable both for e~^e~ and pp collisions, 
but for these first studies we concentrate on the former process type. For each Abelian 
and non-Abelian group we study three different scenarios for the communication between 
the standard sector and the secluded one. We illustrate how to distinguish the various 
characteristics of the models and especially study to what extent the underlying gauge 
structure can be determined experimentally. 
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1. Introduction 

There are basically two ways in which one can envision new physics beyond the standard 
model that can be searched for at future colliders. One possibility is to have theories 
with new heavy particles coupling to the standard model with either the ordinary gauge 
couplings, as in super symmetry, or with couplings of a similar magnitude. This implies 
heavy particle masses, in order to avoid collider constraints. The other possibility, which 
we want to explore in this paper, is that new light particles are ultra-weakly coupled to 
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the standard model particles, because they are not charged under the standard model 
gauge groups. Instead they couple to the ordinary matter through some heavy state which 
carries both SM charges and charges of a new unknown gauge group, also carried by the 
light states. 

There have been several suggestions for theories with this type of secluded sectors 
(sometimes also called hidden valleys or dark sectors), proposing non-conventional new 
physics with unexpected and unexplored signals could show up at current colliders such as 
the Large Hadron Collider (LHC) or a future linear electron-positron collider. 

One example is the so-called hidden valley scenarios by Strassler and collaborators 
1^, ^, ^, ^, 0! where the SM gauge group is extended by a new unspecified gauge 
group G. In the original paper |||] this group is a (/(I)' x SU{N). The new matter sector 
consists of u— particles (where v stands for "valley"), which are charged under the new 
gauge group and neutral under the standard one. The two sectors communicate via higher 
dimensional operators, induced either by heavy particle loops or by a Z' which can couple 
to both sectors. 

An interesting feature of models with secluded sectors is that they naturally give rise to 
dark matter candidates. Likewise, some of the recently proposed dark matter models may 
present hidden sector features. Specific dark matter models developed in the last few years 
such as m, 0, |ll|, 12] suggest the existence of a GeV scale mass dark photon or scalar 
that is introduced to enhance the dark matter annihilation cross section, in order to fit the 



data from PAMELA |lj and originally ATIC although the latter data have later 



been superseded by more precise measurements from FERMI ||T^. Here we will mainly be 
interested in models with dark photons originating from a softly broken ?/(l), which couple 
to standard model particles through so called kinetic mixing with the ordinary photon |17| 
through heavy particle loops in a similar way to the hidden valley scenarios. 

The hidden valley-like theories and the dark matter models mentioned above share two 
features: the enlarging of the standard model symmetries to include a new gauge group G 
and the presence of new light particle sectors that are solely charged under this new gauge 
group. If the new light particles can decay into standard model particles, their existence 
could be inferred from their effect on standard model particle phenomenology. In we 
studied the effects of the new gauge group radiation, specifically the kinematic effects of 
SU{3y radiation from fermions charged under both the SM and the new gauge group on the 
kinematic distributions of visible particles. In this paper we address the issue of discerning 
between different gauge structures. Specifically, we want to outline the differences between 
signatures arising from a secluded sector broken ^(1)' gauge group with a light 7' and 
those arising from a confining SU{N). In both cases we assume there is a mechanism for 
the secluded particles to decay back into the SM. 

In order to distinguish which features of a given model are linked to the gauge group 
structure and which to the other details of the model, we consider different production 
processes and different mechanisms for the decay back into the SM. The various possibilities 
are summarized in Fig. |l]. For the production we consider three mechanisms. In the first 
case the portal to the hidden sector is through kinetic mixing between the SM photon and 
a light ^(1) gauge boson. In the second case we have production via a Z' . In the kinetic 
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mixing of 7-71,, the 7^, is assumed to have a mass around 1-10 GeV and the mixing e is 
assumed to be e ~ 10~^, while in the Z' case, the mass of the Z' would be around 1-6 TeV 
Q. The third case is where the production happens via SM gauge bosons as in |l^. In 
this case the particles are assumed to have both SM charges and secluded sector ones. To 
distinguish them we will call them in the following, to separate them from particles that 
are only charged under the secluded gauge group, which we call q^. We will also assume 
that the particles will decay into a standard model particle / and a secluded sector 
particle g^, i.e. F^ — )• fq^. 

If the particles of the secluded sector are charged under a non-Abelian SU (N) or a 
softly broken Abelian {^(1) with a light 7^,, there will also be additional radiation of gauge 
bosons. In the former case, the 7;-gluons will be connected to the q^s (produced directly or 
via the F^s) and form a confined system which will then hadronize. 

Depending on the nature of the secluded hadrons thus produced, they may then decay 
back into standard model particles through kinetic mixing or a heavy Z' . In both cases 
this decay can be very slow, so much so as to generate displaced vertices and other exotic 
signatures as for example discussed in In the case of 7„ radiation instead, the gauge 
bosons may decay directly back into the SM through kinetic mixing 7f~T; while the qyS 
will not be able to decay back into the SM since they carry the secluded gauge charge. 

Thus, in both the non-Abelian and Abelian cases, we can have models where some of 
the particles produced will decay back into SM particles and some of them will be invisible. 
The questions we want to address is thus how the production of visible particles depends on 
the secluded gauge structure and whether it is possible to tell a non-Abelian and Abelian 
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gauge group apart also when other features of the models are very similar. 

The paper is structured as follows. Section ^ gives a short overview of the general 
model considerations that underlie our studies, with particular emphasis on the produc- 
tion mechanisms that are relevant in various scenarios, and some comments on the decay 
mechanisms that lead to signals in visible distributions. In section |3| we provide a more 
in-depth overview of the new physics aspects that we have implemented in Pythia 8: the 
particle content, the parton showers, notably for the broken U{1) case, the hadronization 
in the secluded sector, and the decay back to the visible one. In Section ^ we proceed 
to describe the phenomenology of the various scenarios, in the context of an e~^e~ lin- 
ear collider. While less interesting than a corresponding LHC phenomenology, it allows 
us to better highlight the relevant features of the model as such. In section ^ we further 
study distributions that could offer a discrimination between an Abelian and a non-Abelian 
scenario for the secluded sector. In Section |^ we summarize our findings, and give an out- 
look. Finally in Appendix |^ we provide information how the simulation of a wide range of 
scenarios can be set-up. 

2. Overview of hidden sector scenarios 

As already mentioned in the introduction there are many different types of models that can 
display hidden sectors and the common feature is that they communicate with the standard 
model through some heavy states. This communication can occur in many different ways 
and we will distinguish three different types in the following: via kinetic mixing, via a 
heavy Z', and via heavy fermions that carry SM charges. 

The common feature of the models we consider is that the SM group SU{3)c x SU (2) l x 
U{1)y is augmented by a new gauge group G. For each scenario we will consider two cases 
- one Abelian where G contains a softly broken ^(1) with a light gauge boson 7^ and one 
non-Abelian where G contains an unbroken SU{N) factor mediated by a then massless gy. 

The particle content consists of particles and/or particles. With particles 
we indicate fermions or scalars (with spin = 1/2, 0, 1) charged solely under the new 
gauge group. With Fy we indicate particles (spin s = 0, 1/2, 1) which may couple to both 
secluded sector and standard model sector. Though in principle one could choose any spin 
assignment among the ones above, we have chosen to analyze the case in which Fy and 
Qy are fermions, except in the case when the qyS are produced from a Fy decay when we 
assume them to be scalars. In addition, in all the scenarios we consider both Fy and qy 
belong to the fundamental representation of the group G. Finally, the G sector charges are 
constrained by anomaly cancellation. For an example see |||]. 

2.1 Kinetic mixing scenarios 

As already alluded to, one way of producing the secluded sector particles is through kinetic 
mixing. In the scenarios we wish to investigate, the SM U(l) (effectively the photon) mixes 
kinetically with a new GeV mass 7' and produces a pair of secluded sector particles, see 
Fig. ||. This mechanism is mostly relevant in the case when the secluded sector contains 
new fermions which are charged only under the new gauge group G. In addition we will 
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in this scenario only consider those cases when the SM particles are not charged under the 
new U(l). Communication between the SM and secluded sectors then only comes from 
kinetic mixing between the standard model U(l) gauge boson and the new gauge boson, 
as described by 

C^n = -\e, {Frf - leFrF^,,. - \e, {F^' • (2-1) 

In the case of two U{1) gauge symmetries {U {l)ixU {1)2) , the non- vanishing mixing e arises 
naturally as one integrates out loops of heavy fermions coupling to both the associated 
gauge bosons [17] so long as there is a mass splitting among them. The relation between 
the size of the mixing and the mass splitting is given by 

where ei and 62 are the gauge couplings of the fermions in the loop to the two U{1) 
gauge bosons, Ai and the new A2 respectively, and and are their masses. In 

general, the U{l)i and the f/(l)2 will not be orthogonal. One may however chose the U{l)i 
generator so that the fermions that are only charged under f7(l)i do not have any charge 
shift, while those that couple to U{1)2 do [17|. 

For the case of non-Abelian groups, Gi x G2 x G3, a mixing can come from the 
spontaneous breaking of the group down to H x f/(l)i x U (1)2- Also in this case the U{l)i 
and the U{1)2 will not be orthogonal, as long as the three couplings associated to the 
unbroken symmetries are different. 

The kinetic mixing mechanism has been used in model that want to describe various 
recent cosmic ray measurements in terms of dark matter models. The most important 
signal here is the positron excess observed by PAMELA |14|. At the same time, any model 
wanting to explain this excess also has to explain the absence of an anti-proton excess 
observed by PAMELA [13| and finally the measurements of the total electron and positron 
flux observed by the Fermi LAT collaboration [|T6[| . The models are set up so that the dark 
matter particles will annihilate into a dark photon or scalar which couples to SM particles 
through kinetic mixing. The mass of the dark matter particle is then determined by the 
scale at which the positron excess is observed, to be of order 0.1-1 TeV. 

In addition, the large positron excess observed also means that there must be some 
enhancement mechanism of the dark matter annihilation cross section. One way to do this 
is to invoke Sommerfeld enhancement^ by introducing a light dark photon or scalar. The 
mass of the dark photon (or scalar) in these models is typically in the GeV range, which 
means that decays into p and 7r° are kinematically suppressed relative to the lepton decays 
and thus also explain the non-observation of any anti-proton excess by PAMELA fistj . 

A recent example of models that fits all these data is given by but there are still 
large uncertainties due to cosmological assumptions such as the dark matter distribution 
and propagation of cosmic particles. 

The dark gauge group Gdark is largely unspecified in these types of models except that 
it must contain a U(l) factor in order for the kinetic mixing with the SM photon. This 



""^Resummation of t-channel exchanges of a new light particle. 
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means that there could also be additional Abelian or non-Abelian factors in Gdark- In 
the following we will consider the cases when Gdark contains an additional U{1), which is 
spontaneously broken giving a massive Z' , or an additional SU{N) factor giving a confining 
force for the secluded sector particles. 

The phenomenology and constraints on these types of models at low energy e~^e~ 



colliders such as Belle, BaBar, DA^>NE, KLOE and CLEO have been studied by pOl, pT 



22, 23 



2.2 Z' mediated scenarios 

The second type of scenarios we want to consider are those that are similar to the original 
hidden valley scenario [Q] with a massive Z' coupling to both SM fermions and secluded 
sector ones. Thus, the processes we are interested in are when SM fermions annihilate into 
the secluded sector Z' which in turn gives a pair of secluded particles, as depicted in Fig. |l]. 

In these types of models it is typically assumed that the Z' acquires a mass by sponta- 
neous symmetry breaking of a U{1) symmetry by a {(p) whereas the origin of the secluded 
sector (/(I) is not discussed. 

The secluded sector particles that the Z' would decay to could be either charged solely 
under the valley gauge group G or charged under G and (parts of) the SM SU{3)c x 
SU{2)l X U{1)y- In the latter case, the particles would on the one hand have to be 
very massive (several hundreds of GeV) due to experimental constraints and on the other 
hand they would be more effectively produced through their SM couplings. Thus we will 
not consider this possibility more here. In contrast the particles charged solely under the 
secluded gauge group could be light with a mass in the 1 — 50 GeV range, thanks to the 
reduced coupling through the heavy Z' . 

As a consequence of the heavy mass of the Z' , the s-channel pair production cross 
section will be peaked at ^fI ~ mz' and be suppressed at an e+e" collider unless ~ Mz'- 
At a hadron collider the production of the Z' would be dominantly on-shell if the overall 
center of mass energy is large enough and there is enough support from the parton density 
functions. 

In the original hidden valley model the secluded sector group also contains a confin- 
ing SU{N). Thus the produced secluded sector particles would have to hadronize into 
hadrons which are neutral under this SU{N). Another possibility is that there is instead 
an additional (/(I) which would instead give radiation of 7's. 

Finally it should be noted that also in this case there is kinetic mixing between the Z 
and the Z', which primarily is important for setting limits on the mass and couplings of 
the Z' from LEP as discussed in Q. 

2.3 SM gauge boson mediated scenarios 

The final type of scenario that we consider are ones where the "communicator" is charged 
under both the SM and new interactions. This scenario and its implementation into pythia 
8 has been described in |jl^ so here we only briefly recapitulate the main features. 

In this model the new heavy communicator particle would be pair produced with 
SM strength, which means that it would have to be quite heavy in order to not have been 
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already seen at colliders. Another consequence is that the communicator would decay into 
a SM and pure hidden sector particle, dubbed qy, so that quantum numbers are conserved. 
In the simple case in which neither q^s nor f-gauge bosons leak back into the SM, as in 
the scenario in [18|, this entails a missing energy signal. 



Also in this case, the secluded sector group can be either Abelian or non-Abelian. In 
both cases we will assume that the produced 7's or hadrons can decay back to SM particles 
through kinetic mixing via loops of the particles or via a Z'. 

2.4 Decays back to the SM 

First of all we mention again the case of secluded particles which are charged both under 
the SM and secluded gauge groups, Fy, which we assume decay according to Fy — )• fqy. 
All other particles produced by either of the mechanisms described above may decay back 
to SM particles as long as they do not carry any charge under the secluded gauge group. 
Essentially these decays will be through kinetic mixing with SM gauge bosons or through 
a heavy Z' as detailed below. 

In the Abelian case, with a light secluded sector 7', the qyS will be stable, but the 
7's that are radiated in connection with the primary hard process will decay back to SM 
particles, 7' — >• //. The strength of the kinetic mixing e, together with the available phase 
space, determines the decay width Ty^j-j. Since the 7' is light, it will mainly mix with the 
standard model photon and thus the branching ratios for different channels will depend on 
the electric charge of the produced SM particles. In essence this means that the decays 
will be similar to a off-shell photon, 7* with the virtuality given by my. We also note that 
if the kinetic mixing is small, the life-time could be so large as to give displaced vertices. 

In the non-Abelian case the secluded sector hadrons may also decay back into the 
SM via kinetic mixing of the 7' with the SM photon or via a heavy Z'. In this case the 
phenomenology will depend on the number of light flavours A^flav in the secluded sector. In 
the following we will assume that Afl^v > 2 and only consider the case when the fundamental 
particles are fermions as in although similar arguments can be made also in the case of 
scalar constituents. Thus, the bound states will be the secluded sector version of mesons, 
baryons and possibly also glueballs. For the decays back to SM particles, it is the meson 
states that are of primary interest and therefore we concentrate on them here. 

With Aflav light flavours, there will be of the order Nq^_^ mesons with a given spin 
out of which approximately Aflav are flavour neutral and can decay back into the SM via 
kinetic mixing or a Z' . The SM decay products will depend on the spin of the secluded 
meson. For a spin zero meson, helicity suppression leads to dominance by the heaviest SM 
particle available whereas for a spin 1 meson it will depend on the couplings to the particle 
mediating the decay, i.e. either to the photon in the case of kinetic mixing or to the Z'. 

The phenomenology will thus depend on the relative production of spin-0 and spin- 
1 mesons and their masses. If the confinement scale A^, in the secluded sector is large 
compared to masses of the lightest secluded sector fermions the situation will be similar 
to QCD. In other words there will be a light spin-0 iTy with mass much smaller than the 
spin-1 py. Thus all p„ will decay to pairs of tt^s and the SM particles produced will be the 
heaviest one available. 
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If At, is of the order of the masses of the hghtest secluded sector fermions then the mass 
sphtting between the spin-0 and spin-1 mesons will be small and thus the spin-1 meson 
will be metastable and instead decay back into the SM, again via either kinetic mixing or 
a Z', but in this case, there not being any helicity suppression, the decay will be similar 
to that of an off-shell photon. Thus in this case there will also be an abundance of leptons 
produced along with hadrons. 

If all constituent masses are much larger than the confinement scale, the lowest lying 
SU{N) neutral states would be glueballs as discussed in We do not discuss their 
phenomenology here. We will also not consider so called quirks which are charged 
both under the SM SU{3)c and a secluded SU{N) with the confinement scale A being 
much smaller than the Fy masses. 

Finally we note that similarly to the Abelian case some of the secluded sector hadrons 
could be metastable and decay back into the detectors with displaced vertices. 

3. Physics in the secluded sector 

For the studies in this article we have developed a framework to simulate the physics of 
a secluded sector. It contains a flexible setup that can be used to study different produc- 
tion mechanisms, perturbative shower evolution scenarios, non-perturbative hadronization 
sequences and decays back into the visible sector. Parts of the framework were already in 
use for our previous study ||l8| but significant new capabilities have been added. These are 
available starting with Pythia 8.150. The physics content will be described in the follow- 
ing, while technical details on how to set up a variety of scenarios is outlined in Appendix 
A. The studies presented in this article only give a glimpse of the possibilities. 

3.1 Particles and their properties 

The key aspect of a scenario is that of the valley gauge group G, which we allow to be 
either U{1) or SU (N). The gauge bosons of these groups are named 7i, and g^, respectively. 
The former can be broken or unbroken, i.e. 7^ can have a mass, while the latter is always 
unbroken so that remains massless. 

The rest of the particles, i.e. the "matter" content, fall into two main categories: those 
charged under both the SM and the v sector, and those that are pure f-sector particles. 

For the doubly charged ones, dubbed F^, 12 particles are introduced to mirror the 
Standard Model flavour structure, see Tab. |2| in the appendix. Each F„ particles couples 
flavour-diagonally to the corresponding SM particle. In addition to its SM charges, it is 
also put in the fundamental representation of G. For U{1) the charge is taken to be unity, 
while for SU{N) the "charge" is Cf = (A^^ — 1)/{2N) while pair production cross sections 
obtain a factor of enhancement. Although the name suggests that the Fy are fermions, 
they can be spin 0, 1/2 or 1 particles. If the Fy particles have spin 1 then their production 
cross section depends also on the presence or not of an anomalous magnetic dipole moment. 

The valley secluded sector further contains a purely G interacting sector. At the parton 
level this consists of q^s, belonging to the fundamental representation of G. The name is 
introduced to reflect the similarities with the quark in QCD. The qy particle is stable and 
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invisible to SM interactions. Its spin, or 1/2, is adapted to the choice of spin made for 
Fy, in case the scenario allows for Fy — )• fq^ decay, where / is a SM particle. The spin 
structure of the Fy — )• fq^ decay is currently not specified, so the decay is isotropic. 

In the G = U{1) scenarios only one q^ is assumed to exist. Fy decays, if allowed 
kinematically, are fiavour diagonal, Fy — )• f^qy, with a common (Yukawa) coupling strength. 
Given that both the FyS and the q^s have a unit of U{1) charge, they can radiate 7^ 
gauge bosons. If U{1) is unbroken the 'jy is massless and stable. For a broken symmetry, 
G =^(1), the 'jy can decay back to a SM fermion pair through the mechanisms discussed 
in the previous Section ^. For kinetic mixing or decay via a Z', branching ratios by default 
are assumed to be proportional to the respective fermion coupling to the photon, whenever 
the production channel is allowed by kinematics. The 'jy decay can be either prompt or 
displaced. 

If instead G = SU{N), the massless Qy gauge bosons are self-interacting, such that the 
parton shower will also have to allow for Qy — )• gygy splittings, with no equivalence in the 
U{1) case. The self-interactions also lead to confinement, like in QCD. In Section |3.3| below 
we will explain how the resulting picture can be described in terms of "strings" stretched 
from a qy end via a number of intermediate ^^.s to a end. The string can break, by the 
production of new qyqy pairs, to produce a set of u-mesons formed by the qy of one break 
and the qy from an adjacent one. To first approximation these i;-mesons would be stable, 
and so the whole u-hadronization process would be invisible. One would not even have the 
kind of indirect recoil effects that the u-shower can give. If kinetic mixing or decay via a 
Z' is assumed, it would again be possible to let the u-mesons decay back to a SM fermion 
pair. 

With only one qy species there would only be one kind of u-mesons, and so the choice 
would be between two extremes: either all the energy deposited in the hidden sector decays 
back to be visible, or none of it. The more interesting scenarios — e.g. in terms of offering a 
bigger challenge to sort out what is going on — are the ones where only part of the f-mesons 
can decay back. Therefore a variable number A^flav of separate qy fiavours are assumed to 
exist (at most 8 in the current implementation). This gives A^^^^^^ different possible f- meson 
flavour combinations, out of which only iVfjav are flavour-diagonal and thus able to decay 
back into the SM sector. It would be possible to assign individual masses to the qyS and 
u-mesons, but for now we assume one common qy "constituent" mass and one common 
i;-meson mass, twice as large as the former. 

By analogy with QCD two separate spin states are assumed, denoted TTy and py. For 
now mass splitting is taken to be small, such that py — >• TTyiTy is kinematically forbidden, as 
is the case in QCD for the s and heavier quarks. The decay of the flavour-diagonal mesons 
is different in the two cases: by helicity (non) conservation the TTy couplings to a pair of SM 
fermions / provides an extra factor mj , an addition to the squared charge and phase space 
factors factors present for the py mesons. 

In the conflning SU (N) case also a I'-glueball is introduced. It is only rarely used, 
to handle cases where the invariant mass of the invisible-sector fragmenting system is too 
large to produce one single on-shell f-meson and too small to give two of them. Then it is 
assumed that an excited v-meson state is produced, that can de-excite by the emission of 
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these invisible and stable gyg^ bound states. 

In summary, by default the v-particles with no SM couplings are not visible. Their 
presence can only be deduced by the observation of missing (transverse) momentum in the 
event as a whole. On top of this we allow two different mechanisms by which activity can 
leak back from the hidden sector. The first is the Fy — )• fq^ decay and showers from the 
Fy and qy, in the scenario in which Fy has both SM charges and G charges, as discussed in 



our previous article [18|. The second is the decay of SM gauge bosons produced through 
mixing by the G group gauge bosons in the kinetic mixing case, either the massive ^y for 
1/(1) or the diagonal t;-mesons for SU{N). 

3.2 Valley parton showers 

Parton showers (PS) offer a convenient approximation to higher-order matrix elements, 
which by the use of Sudakov form factors contain a resummation of virtual corrections 
to match the real emissions [26|. For the current studies, the Pythia p_|_-ordered parton 



showers are extended to the secluded sector, and the approach used to take into account 
massive radiating particles |28| must, for the ?/(!) scenario, be extended to the case where 



also the radiated gauge boson is massive. This section gives a summary of the showering 
framework, with emphasis on aspects new to this study (relative to ]18[|). 

In the most general case, final-state QCD, QED and valley radiation are interleaved 
in one common sequence of decreasing emission p± scales. That is, emissions of a SM g/j 
or a hidden 'jy/gy can alternate in the evolution of a Of course any of the related 
charges can be zero in a specific process, in which case the following expressions simplify 
accordingly. For the i'th emission, the p± evolution starts from the maximum scale given 
by the previous emission, with an overall starting scale p±Q set by the scale of the hard 
process, or of the decay in which the radiating particle was produced. Thus the probability 
to pick a given p± takes the form 

dV _ fdVQCB ^ dpQED ^ dPsecA /_ f"^^-' / dpQCD ^ dT^QED ^ dPsecA 



dp± V dp_i_ dp± J \ Jp^ \ dp'_L dp'j_ dp'j 

(3.1) 

where the exponential corresponds to the Sudakov form factor. Implicitly one must also 
sum over all partons that can radiate. 

To be more precise, radiation is based on a dipole picture, where it is a pair of partons 
that collectively radiates a new parton. The dipole assignment is worked out in the limit 
of infinitely many (hidden or ordinary) colours, so that only planar colour flows need be 
considered. 

Technically the total radiation of the dipole is split into two ends, where one end acts as 
radiator and the other as recoiler. The recoiler ensures that total energy and momentum 
is conserved during the emission, with partons on the mass shell before and after the 
emission. Each radiation kind defines its set of dipoles. To take an example, consider 
qq —7- FyFy, which proceeds via an intermediate s-channel gluon. Since this gluon carries 
no QED or hidden charge it follows that the FyFy pair forms a dipole with respect to these 
two emission kinds. The gluon does carry QCD octet charge, however, so FyFy do not 
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form a QCD dipole. Instead each of them is attached to another parton, either the beam 
remnant that carries the corresponding anticolour or some other parton emitted as part of 
the initial-state shower. This means that QCD radiation can change the invariant mass of 
the FyFy system, while QED and hidden radiation could not. When a 7 or 7^ is emitted 
the dipole assignments are not modified, since these bosons do not carry away any charge. 
A g or gy would, and so a new dipole would be formed. For QCD the dipole between Fy 
and one beam remnant, say, would be split into one between the Fy and the g, and one 
further from the g to the remnant. For the secluded sector the FyFy dipole would be split 
into two, Fygy and gyFy. As the shower evolves, the three different kinds of dipoles will 
diverge further. 

Note that, in the full event-generation machinery, the final-state radiation considered 
here is also interleaved in p± with the initial-state showers and with multiple parton-parton 
interactions [p9|]. 

If the Fy fermion is allowed to decay into a SM and a hidden particle, one must also 
consider the hidden radiation from the hidden particle. 

There is a clean separation between radiation in the production stage of the FyFy pair 
and in their respective decay. Strictly speaking this would only be valid when the Fy width 
is small, but that is the case that interests us here. 

In the decay Fy — >• fqy the QCD and QED charges go with the / and the valley one 
with Qy. For all three interactions the dipole is formed between the / and the qy, so that 
radiation preserves the Fy system mass, but in each case only the relevant dipole end is 
allowed to radiate the kind of gauge bosons that goes with its charge. (Strictly speaking 
dipoles are stretched between the / or qy and the "hole" left behind by the decaying Fy. 
The situation is closely analogous to t — )• bW~^ decays.) 

The number of parameters of the hidden shower depends upon the scenario. In the 
case of the interleaved shower, there are only two, the most important on being one the 
coupling strength Oy, i.e. the equivalent of Ug. This coupling is taken to be a constant, i.e. 
no running is included. 

From a practical point of view it is doubtful that such a running could be pinned down 
anyway, and from a theory point of view it means we do not have to specify the full flavour 
structure of the hidden sector. The second parameter is the lower cutoff scale for shower 
evolution, by default chosen the same as for the QCD shower, p^min — 

0.4 GeV. 

3.2.1 Shower kinematics with massive hidden photons 

Showers are expected to reproduce the soft and collinear behaviour of (leading-order) ma- 
trix elements (MEs), but there is no guarantee how trustworthy they are for hard wide- 



angle emissions. Therefore various correction techniques have been developed [26|. The 
technique we will use here is to generate trial emissions according to the PS, but then use 
the weights ratio ME/PS to accept emissions, i.e. PS times ME/PS equals ME. For this 
re-weighting recipe to work, obviously the ME weight has to be below the PS one, but the 
difference should not be too big or else the efficiency will suffer. It should also be noted that 
the ME/PS ratio is evaluated without including the Sudakov form factor of the shower, 
while the shower evolution itself does build up the Sudakov. By the veto algorithm it then 
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follows that the ME expression is exponentiated to provide the kernel of the Sudakov [pO|, 



a technique nowadays used as a key ingredient of the POWHEG approach |31]. The choice 
of shower evolution variable lives on in the integration range of the Sudakov, but for the 
rest the PS expressions disappear in the final answer. 

In the past, this approach has only been developed for the emission of a massless gluon, 
however, and we now need to generalize that to an arbitrary combination of masses. A 
technical task is to recast the ME and PS expressions to use the same phase space variables, 
such that the ratio is well-defined. 

We follow the existing approach of mapping the PS variables onto the ME ones. Below 
we therefore introduce the ME three-body phase space, subsequently how the PS variables 
populate this phase space, and finally how the presence of two shower histories can be 
taken into account. 

Consider a dipole of invariant mass niQ, consisting of two endpoint partons 1 and 2, 
with nominal masses mi and m2- Assume that a shower emission occurs from the parton- 
1 dipole end, generating a new particle 3 with mass m^. This implies that there was 
an intermediate off-shell state 13 with mass m-13. That is, the kinematics to describe is 
Po ~^ Pi3 + P2 Pi + P?,+ P2- Averaging over the angular orientation of events, the MEs 
can be written in terms of the Xi = 2pipo/m1 and the ri = mf/ml variables, where the Xi 
reduce to energy fractions in the dipole rest frame, with normalization xi + X2 + x^ = 2. 
This means there are only two free independent variables, traditionally xi and X2. 

The PS is instead described in terms of the Pj_cvoi ^"^^ ^ variables. In the soft and 
collinear emission limit these are well defined, but away from these limits different possi- 
bilities could be contemplated. Our choice is such that 

ml, = ml + , (3.2) 

or 

Pievol = z{l - z){ml^ - ml) . (3.3) 
By standard two-body kinematics for pQ — )• pi, + p2 it follows that 

X2 = 1 = 1 + ^2 - ri3 , (3.4) 

and thus xi + x-, = 2 — X2 = 1 + ?'i3 — T2. If mi = m, = one would further require that 
z = xi/{xi + X3). Taken together, this is enough to specify the three four-vectors p2, p^i^ 
and ^'3'^^ up to three angles. These are chosen at follows: in the po rest frame parton 2 
is assumed to keep its direction of motion when mi — )• mi3, while 1 and 3 are selected to 
have an flat distribution in the azimuthal angle around the 13 direction, which is parallel 
with the 1 direction before the emission. 

The kinematics for the case with massive partons 1 and 3 can then be constructed 
from the massless four-vectors as 

pi = {1 - h)pf^ + kjp (3.5) 
P3 = (1 - h)pP + kipf (3.6) 
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mjg - Ai3 ± {ml - ml) 
2mi 



Ai3 = Y {ml^ — ml — 7713)2 — Amlm"^ (3.8) 

The physics content is that the directions of partons 1 and 3 in the P13 rest frame are 
retained, while their three-momenta are scaled down by a common factor sufficient to put 
the two partons on their mass shells. Since mis is not changed by the operation it is 
necessary that mi3 > mi + ms for the rescaling to work. 
The rescalings imply that 

= = (1 - ki)z + A;3(l -z) = il-ki- h)z + ks , (3.9) 

Xl +X3 2 - X2 

and thus 

^ ^ ks) ■ (3.10) 



I - ki - ks \2 - X2 

Now we need to find the Jacobian to translate the shower emission rate from the 
(Pj_evoi' ^) space to the (xi, X2) one. Note that ml^ = mo(l — X2) + m| is independent of xi, 
and thus so are ki and k^. Therefore only the "diagonal" terms dp']_^^^Jdx2 and dz/dxi 
are needed. 

The shower emission rate itself is 



dpjevoi 2dz 



(3.11) 



Here an overall coupling factor Cp 01^/211 is omitted for simplicity. Also the Sudakov form 
factor is omitted, as already motivated. The z-dependent part may seem unfamiliar, but 
is an upper approximation to the more familiar q ^ qg splitting kernel (1 + z'^)/{l — z), 
where the difference between the two is absorbed into the ME/PS weighting. 
Put together, the shower emission rate translates into 



dpievoi _ d(mf3 -m'f) _ dxa 



^'ievoi "^13 -^1 1 - X2 + ra - ri 
2dz 2 dxi 



(3.12) 



l-z (1 - /ci - /C3)(2 - X2) 1 1 

2dxi 



X3 - ki{xi +X3) 

When ms — t- 0, and hence ki — )■ 0, this simplifies to the familiar expression |28] 



(3.13) 



Wps, = ^1^= '^^'^-^ ^ . (3.14) 

^'levol 1 ~ ^ (1 - X2 + r2 - ri)X3 



If only parton 1 can radiate, as in — )■ + / — g^, + 7^, + /, we are done. The 
fact that the MEs also contain a contribution from 7„ emission off the does not change 
the picture, since that does not introduce any new singularities, and empirically the PS 
expression provides a valid upper limit. 
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For the radiation F^F^ — t- F^F^j^ the sum of the two possible shower emissions are 
needed to match to the fuU MEs. Alternatively, and more conveniently, the ME expression 
can be split into two parts, each to be compared with only one shower history. This split 
is done in proportion to the respective propagator, i.e. assumed emission off parton i is 
proportional to l/{mf^ — m^). The relative probability for parton 1 to radiate thus is 

_ m^g -ml _l-xi + ri-r2 



(mfg - mf ) + (m|3 - m|) X3 
so that the ME weight to be associated with this dipole end is 



(3.15) 



Wme.i = Pi— , dxi dx2 . (3.16) 

(To dxi dX2 

Thus we arrive at the ME/PS correction factor 

„ Wme,i (1 - xi +ri - r2)(l - a;2 + r2 - ri) 1 da 
Jti = — 



Wps,i 2 (To dxi dx2 

X3 - ki{xi +X3) 



(3.17) 



X3 

All the explicit dependence on 771,3 is located in ki in the last term, but obviously implicitly 
the whole kinematics setup is affected by the value of 7713. 

The matrix elements for the radiation off F^F^ are calculated with them as stable 
final-state particles. This means that, to preserve gauge invariance, they must be assigned 
the same mass. On the other hand, since they are supposed to decay, we allow them to have 
a Breit-Wigner mass distribution. To resolve this discrepancy, the real kinematics with two 
different masses is shifted to a fictitious one where F^ and F-u have the same mass, and it is 
this fictitious one that is used in the three-parton matrix-element evaluation. As a guiding 
principle, the F^ and F^ three-momenta are kept unchanged in the F^Fy rest frame, and 
only energy is shuffled so as to equalize the masses. Denoting the average mass m, the 
conservation of three-momentum implies that 



"^12 —2 / ("^12 ~ "^1 ~ fn^Y — 4777,^777,2 

4?7T,^ 



'""12 



which gives 

= -^^^-^ - . (3.19) 



777^ + 7772 ("^1 ^ "'2)'^ 

~2 4m? 



''12 

As above, the modified four-vectors Pi and P2 can be written as linear combinations of the 
original ones, with the constraints p\=P2~ giving the solution 

2 2 

Pi=Pi+ "'l~2^^ (Pi + P2) , (3.20) 

2 2 

P2=P2~ I 2 iPi + P2) . (3.21) 

/777^2 

This translates into identical relationships for the modified matrix-element variables xi 
and X2 in terms of the original xi and X2 ones. 
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Figure 2: The Feynman diagrams for the production. 



3.2.2 Matrix element for radiation in production 

The implementation of the has required the calculation of matrix element corrections 
IMI^^ ^ for the pair production process // — )• F^F^'yv described in Fig. |2l This has 
required the generalization of the matrix element corrections in [^] to the case of a massive 
photon: 

= (l-4ri)^/2 , (3.22) 

^ p fx 

ff^KF.y. - (rs + 2n)(-l + An) ^ + 



1^1 "V(l-Xl)^ (1-^2)^ 

— 1 + Sri — xo — 1 + Sri — xi 
H H 

1 — Xi 1 — X2 

^^Jl^^Il±^d±±pl + 2. (3.23) 

(1 - Xi)(l - X2) 

Here ri = r2 = /rn^ and = rn^^^/rn^. (Expressions for ri ^ r2 have also been obtained 
but, by the preceding trick, are not needed.) Coupling constants have been omitted, as 
discussed before for the shower. Furthermore, to simplify calculations, the process is taken 
to proceed via the exchange of a scalar particle instead of a spin 1 gauge boson. The 
\M\'^j.j^p p ^ spin information, relevant for decay angular distributions, will be lost this 



way. Effects are known to be minor for the ME correction ratio [28|. As an illustration, 
the above expression reduces to {x\ + ^^/((l — xi){l — X2)) + 2 for ri = r2 = = 0, 
where the first term is the familiar expression for e'^e" — )• 7*/Z* — )• qq, and the second 
finite term comes in addition for a spin exchanged particle. 

3.2.3 Matrix element for radiation in decay 

The matrix elements corresponding to — ?• Qv/jv are 



M||^_ . = (1 - ri + r2 + 2Q2)V(l-r-i-r2)2-4rir2 (3.24) 



|M 



F,^g„/7„ _ (?'3/2 + 2rf + r2r3/2 + ^2^3 - 2ri - rir3/2 - 2rir2 - 4rig2) 



|M|2 (i + ^2-ri-X2)2 



+ 
+ 



-2 + 2r| + 2rf + 2r2r3 - 4^2 + 2g2?'3 + 4q'2?^2 - 4rir2 - 4rig2) 

(1 + r2 - ri - X2)(r3 - X3) 
-2-r3/2-2r2-4g2 + 2ri) 
(1 + r2 - ri - X2) 
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Figure 3: The Feynman diagrams for the Fy — > qv qjl decay. 



(-2 - rs - 2r2 - rirj, - 4^2 - 2g2'^3 + 2ri + t\t-^ 



[rs - xs) 

where ri = m'^^/mp, r2 = m'j/mp, = m^^/m^p and q2 = mf/rriF = The calcula- 

tion has been done for the specific choice of and / being fermions, and a scalar, but 
again the result should be representative also for other spin choices. 

3.3 Hidden sector hadronization 

If the G group is the unbroken SU{N), the gauge boson is massless and the partons are 
confined. The picture therefore is closely similar to that of QCD, and we will use exactly 
the same framework to describe hadronization physics as in QCD: the Lund string model 

For the hidden sector, the model is most easily illustrated for the production of a 
back-to-back q^q^ pair, with the perturbative emission of additional g^s neglected for now. 
In that the partons move apart, the physical picture is that of a ii-colour flux tube 

being stretched between the Qy and the Qy. If the tube is assumed to be uniform along its 
length, this automatically leads to a confinement picture with a linearly rising potential, 
V{r) = nr. 

In order to obtain a Lorentz covariant and causal description of the energy flow due 
to this linear confinement, the most straightforward approach is to use the dynamics of 
the massless relativistic string with no transverse degrees of freedom. The mathematical, 
one-dimensional string can be thought of as parameterizing the position of the axis of a 
cylindrically symmetric flux tube. 

Now consider the simple qyqy two-parton event further. As the qy and qv move apart 
from the creation vertex, the potential energy stored in the string increases, and the string 
may break by the production of a new q'^q'^ pair, so that the system splits into two colour 
singlet systems qyq'y and q'^qy. If the invariant mass of either of these systems is large 
enough, further breaks may occur, and so on until only u-mesons remain. A system with n 
primary mesons thus requires n — 1 breaks qv,iqv,i to produce a chain of mesons qyqy^i, 
Qv,i(lv,2, qv,2qv,3, ■ ■ ■ , Qv,n-iQv Stretching from the q^ end to the ijy one. 
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The flavour of each qv,iqv,i is supposed to be a random choice among the A'flav different 
flavours. Since all are taken to have the same mass, for now, they are also produced at 
the same rate. This thus also goes for the N^.^^ different u-meson flavour combinations. If 
the are fermions then both pseudoscalar and vector t;-mesons can be produced, vr^ and 
Pv. Again disregarding possible effects of a mass splitting, simple spin counting predicts a 
relative production rate 7r„ : p„ = 1 : 3. 

The possibility of higher excited states is disregarded, as is known to offer a good 
approximation for the QCD case. Also f-baryon production is left out, which is a 10% 
effect in QCD. For a generic SU{N) group a u-baryon needs to consist of v-quarks. This 
should lead to exceedingly tiny rates for A^ > 3, while N = 2 could offer a more robust 
u-baryon production rate. 

The space-time picture of the string motion can be mapped onto a corresponding 
energy-momentum picture by noting that the constant string tension implies that the v- 
quarks lose a constant amount of energy per distance traveled. The different breaks are 
space-like separated, but two adjacent breaks are constrained by the fact that the string 
piece created by them has to be on the mass shell for the ?;-meson being produced. The 
space-like separation implies that the fragmentation process can be traced in any order, 
e.g. from one of the endpoints inwards, while the constraint implies that there is only one 
kinematical degree of freedom for each new iJ-meson. Typically it is chosen to be z, the 
light-cone momentum fraction that the new iJ-meson takes from whatever is left in the 
system after previously produced u-meson have been subtracted off. 

By symmetry arguments one arrives at the Lund-Bowler shape of the z probability 
distribution [33| 

m^-^il-zre.p[-^^ , (3.26) 

where rrimy ~ 2mg^ is the mass of the produced w-meson. The equation contains two free 
parameters, a and b. Roughly speaking, these regulate the average rapidity spacing of the 
v-mesons, and the size of the fluctuations around this average. While a is dimensionless, 
b is not, which means that it becomes necessary to adjust b as niq^ is changed. For 
instance, assume that the qy mass is related to the strong-interaction scale A„. Then, if 

iTiq^ , TTirrii, and the collision energy are scaled up by a common factor, we would want to 
retain the same rapidity distribution of produced mesons. This is achieved by rewriting 
6m^^ = (&"iq^ ) (m^^ /m^^ ) = 6'(m^^/m^^), where now b' can be assumed constant. 

In additional to fluctuations in the longitudinal fragmentation, it is assumed that each 
new q'^q[j pair produced when the string breaks also carries an opposite and compensating 
transverse momentum component. The p± of the qv,i-iqv,i meson is then given by the 
vector sum of its two constituent p± values. The pair p± naturally arises in a tunneling 
production process, which also leads to a Gaussian p± distribution. The width a of this 
Gaussian again should scale like A^,, so we rewrite as cr = {c^/mq^)mq^ = a'mq^. When the 
i;-mesons are allowed to acquire a p_\_ it should be noted that the in eq. ( |3.26|) must 
be replaced by m\^^ = m^^ + p\. 

In lack of further knowledge, it is convenient to assign b' and a' values by analogy with 
standard QCD. To be more specific, we have in mind something like the s quark, with a 
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bare mass of the same order as A. For heavy quarks, hke c and b in QCD, tunnehng is 
suppressed, and the framework would have to be further modified. To assess uncertainties 
in a scenario, it would make sense to vary b' and a' values over some range, say a factor of 
two in either direction. 

So far, the emission of ^^.s has been neglected. When it is included, more complicated 
string topologies can arise. Like in QCD, the complexity is reduced by using the planar or 
large- iVc limit |34]. In it a w-gluon is assigned an incoherent sum of a (f-)colour charge and 



a different anticolour one. In a branching qy — QyQy the initial colour is taken away by 
the Qy and a new colour-anticolour pair is stretched between the final qy and gy. Similarly 
Ov — ^ dvdv is associated with the creation of a new colour. That way partons nearby in 
the shower evolution also come to be colour-connected. This leads to a picture of a single 
string, consisting of several separate string pieces, stretching from one qy end to the gy it 
shares one colour with, on to the next colour-related gy, and so on until the qy string end 
is reached. Several separate string pieces could have formed, had perturbative branchings 
dv QvQv been included, but, as in QCD, gy — )• qyqy should be rare both in relation to the 
more singular gy — t- gygy and in absolute terms. 

The motion of a string with several gluon kinks can be quite complicated, but it 
is possible to extend the fragmentation framework of a single straight string also to the 



more complex topologies [35|. Basically the string will break up along its length by the 
production of new q'yq'y pairs, with two adjacent breaks correlated in such a way that the 
u-meson produced between them is on the mass shell. Sometimes the two breaks will be 
on either side of a gy string corner. 

One of the key virtues of the string fragmentation approach is that it is collinear and 
infrared safe. That is, the emission of a gluon disturbs the overall string motion and 
fragmentation vanishingly little in the small- angle/energy limit. Therefore the choice of 
lower cut-off scale for parton showers is not crucial: letting the shower evolve to smaller and 
smaller scales just adds smaller and smaller wrinkles on the string, which still maintains 
the same overall shape. 

The complete f-string fragmentation scenario contains a set of further technical details 
that are not described here. The key point, however, is that essentially all of the concepts of 
normal string fragmentation framework can be taken over unchanged. The one new aspect 
is what to do when the invariant mass of the hidden-valley system is too large to produce 
one single on-shell u-meson and too small to give two of them. As already explained, then 
the emission of f-glueballs is used to balance energy-momentum. 

3.4 Decays back into the SM sector 

Disregarding the trivial direct decay Fy — t- fqy, the main decay modes back into the SM 
are through jy kinetic mixing or Z' decay. For G =J/{1) the jy therefore can decay to SM 
particles with the same branching ratios as a photon of corresponding mass, i.e. oc ej^Nc, 
with A'^c = 1 for leptons. For G = SU{N) only the flavour-diagonal mesons can decay, 
either with a 7^ or a Z'. (The former would imply that G = SU{N) x U{1), which would 
require some further extensions relative to the scenarios studied here.) A py meson, with 
spin 1, could have the same branching ratios as above, or slightly modifled depending on 
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the Z' couplings. A vr^, meson, with spin 0, would acquire an extra helicity factor that 
would favor the heaviest fermions kinematically allowed. Should the vr^, be scalar rather 
than pseudoscalar there would also be a further threshold suppression, in addition to the 
phase space one. 

The decay back into the standard model would be accompanied by normal QED and 
QCD radiation, where relevant. Quarks and gluons would further hadronize, as described 
by the normal Lund string model. That model is not carefully set up to handle different 
exclusive states if the 7„ or Pv/t^v mass is very low, of the order 1 or 2 GeV, but should 
be good enough as a starting point. For studies that zoom in on one specific mass, more 
carefully constructed decay tables could be used instead. 

4. Analysis of the different scenarios 

The tools described above allow us to simulate several different setups. We concentrate 
on the phenomenology of the six scenarios listed in Table |l|. Three different production 
mechanisms are involved: s-channel pair production via kinetic mixing with the light 7„ 
(KM^^), s-channel pair production mediated by a Z' (M^/) and s-channel pair production 
via SM gauge bosons (SM) and the Fy particles. For each of these production mechanisms 
an Abelian setup and a non-Abelian one are considered, labeled by A and NA respectively. 
Note that the Abelian/non- Abelian group we refer to in the following analyses correspond 
to the G gauge group, not to the production mechanisms. In the Abelian case G =^(1), 
while in the non-Abelian case G = SU{3). 
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Table 1: The six scenarios studied. 

The phenomenology of the six scenarios is a function of the pair production cross 
section, which will in general depend upon the specific model realization of each setup. 
In particular, for the KM scenarios, the cross section will depend upon the size of the 
kinetic mixing parameter e, while for the Z' mediate ones on the mass of the Z' and on its 
couplings to the SM particles and to the u-quarks. The analysis is performed on per-event 
distributions, so as to factor out this model dependence. Assuming the same number of 
events are produced, the phenomenology of the setups will depend upon a different number 
of parameters. For the KMA^^ and the AMz' one must fix the Qy masses, the jy mass and 
the coupling constant ay, while for the SMA production one must also fix the Fy 

masses. In the corresponding KMNA^„ , AMz' and SMNA one must fix the meson masses, 
but these will be connected to the qy masses chosen, and furthermore the gy remains 
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massless. We select a scenario in which rriq^ ~ A, so that 171-,^^ ~ i^pv- This in turn ensures 
(as already described in section |2.4| ) that meson decay into SM lepton is not supressed. 
For simplicity, in the following analysis we assume only one mass for all tJ-quark flavours, 
and only one common TTy/p^ mass m^^jp^ = 2mq^. One additional simplification in the 
following analysis, is that for the SM cases we assume the pair production of one single 
belonging to the standard model doublet with no consideration for anomaly cancellation 
issues. In the non-Abelian case we have assumed simple proportions 1 : 3 for iTy : 
production from fragmentation, which comes from spin counting when the has spin 1/2. 
The branching ratios of the decays to standard model particles are fixed by the kinetic 
mixing mechanism. 

We concentrate on the phenomenology of the six setups at an e~^e~ collider with 
center-of-mass (CM) energy of 800 GeV. A similar study for pp colliders like the LHC is 
also possible, and obviously more relevant in the near future, but makes it less transparent 
to compare and understand the properties of the models. Bremsstrahlung corrections have 
been included, and we shall see that these can give a non- negligible effect, whereas the 
machine-specific beamstrahlung has not. All of the figures in this section are based on a 
Monte Carlo statistics of 10000 events. 

As a consequence of the e~^e~ collider choice, the events have a spherical symmetry 
rather than a cylindrical one, i.e. are described in terms of particle energy and {9, (p) 
variables rather than in terms of Et and (77, </>). The jet clustering algorithms are thus 
determined by the spherical topology and we primarily use the Pythia built-in Cluster Jet 
Jade algorithm The Jade algorithm is geared towards clustering objects nearby 

in mass, and so for clustering a variable number of fixed-mass Jv/t^v/Pv systems it is more 
relevant than clustering e.g. in transverse momenta or angles. 

4.1 Basic distributions 

The phenomenology of the six different setups is discussed in detail in dedicated subsections. 
We here would like to discuss the general features of the secluded sector signals and to 
introduce the observables we focus on. 

One of the benefits of a Monte Carlo simulation is that one may look behind the scene, 
to study also the distributions of the invisible secluded sector particles. These can then 
be compared with the SM particle distributions to determine which features are governed 
by the the secluded sector dynamics, and which come from the decays to the SM. In 
this spirit. Fig. Q shows the number of 7^, gauge bosons emitted per event in the AM^/ 
case. On the left we highlight the ay dependence, on the right the m^^ dependence. Not 
unexpectedly, the number of 7^, increases almost linearly with a^, up to saturation effects 
from energy-momentum conservation. 

Compare this distribution with the corresponding non-Abelian AM^/ case in Fig. |5| for 
the flavour diagonal t^v/Pv Again the number of grows with q„ ^, but the number of 
t;-mesons does not primarily reflect this dependence. Instead the number of u-mesons 

^Note that, in this case, the emission rate g„ —>■ qvQv is proportional to Cpdv, with Cf = {N^ — 1)/{2N), 
and Qv — > QvQv to Na^. 
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Figure 4: AMz'- the number of gauge bosons emitted per event. On tlie left we Irigliliglit tlie 
ay dependence, while on the right the m^^ dependence. On the left side ruq^ — 50 GeV, m-^^ = 10 
GeV and = 0.1, 0.4, 0.6, while on the right side m~^^ = 2, 6, 10, 20 GeV and the coupling is fixed 
at a„ — 0.4. 




Figure 5: NAMz''- the number of flavour diagonal TTy/py mesons emitted per event. On the left 
we emphasize the ay dependence, while on the right the 'rij^^/p^ dependence. On the left the meson 
mass is fixed at m^^/p^ — 10 GeV and the coupling varies among ay — 0.1,0.4,0.6, while on the 
right side m^^/p^ = 2,6,10,20 GeV (which in turn implies ruq^, = 1,3,5,10 GeV), the coupling is 
fixed at a„ = 0.4 and the number of flavours is A^flav — 4. 



produced by string fragmentation rather reflects the masses of the T;-quarks (and thereby 
of the mesons) and the fragmentation parameters, see Fig. |5| right plot. Specifically, even 
with set to zero for the perturbative evolution, there would still be non-perturbative 
production of f-mesons from the single string piece stretched directly from the to the 
Qy. With a-u nonzero the string is stretched via a number of intermediate gluons that 
form transverse kinks along the string, and this gives a larger multiplicity during the 
hadronization. 

Comparing the number of 7^, in Fig. ^ with the corresponding distributions in the other 
Abelian setups, KMA^^ and SMA, in Fig. 12 and Fig. 14 respectively, the two KMA^^ and 
AM^/ setups produce similar distributions, while the SMA produces much fewer 7^,. The 
SMA difference is due to the more complicated kinematics, where the electrons from the 
— >• eq^ decays take away a large fraction of energy and momentum that then cannot be 
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Figure 6: AMz' , NAMz''- dependence of the overall number of charged particles emitted per 
event in the Abelian (left) and non-Abelian (right) case. For the Abelian cases ruq^ — 50 GeV and 
m^^ — 10 GeV, for the non-Abelian cases ruq^ = 5, vHt^^i — 10 GeV. 

used for 7^ emissions. 

The average charged multiplicity of an event, Fig. ^, will be directly proportional to 
the number of ^v/t^v/ Pv produced. The trends from above are thus reproduced, that the 
non-Abelian multiplicity varies only mildly with a^, while the variation is more pronounced 
in the Abelian case. The constant of proportionality depends on the ^v/t^v/ Pv mass, with 
more massive states obviously producing more charged particles per state. This offsets 
the corresponding reduction in production rate of more massive ^v/t^v/ Pv, other conditions 
being the same. Similarly the number of jets should be proportional to the number of 



Iv/t^v/ Pv emitted, see further Sec. [O 



Without an understanding of the ^v/t^v/ Pv mass spectra, the mix of effects would make 
an determination nontrivial, especially in the non-Abelian case. Even with a mass fixed, 
e.g. by a peak in the lepton pair mass spectrum, other model parameters will enter the 
game. One such parameter is the number N^^^ of flavours. Since only l/A'^flav of the 
T^v/ Pv would decay back into the SM the visible energy is reduced accordingly. With all 
having the same mass, the relation (-E'visibie)/-E'cm = l/-^flav works fine to determine N^^,^, 
but deviations should be expected for a more sophisticated mass spectrum . Furthermore 
the 7r„ : p„ mix, with different branching ratios for the two, needs to be considered. If 
the 7r„ fraction is large, the number of heavy leptons and hadrons produced may increase 
substantially, see 

In Fig. we show the energy spectra of the hidden sector 7^, and Pv/t^v Note the 
difference between the NAM^/ setup and the KMNA7„ one. This is due to the difference 



in the amount of initial-state radiation in the two cases, as discussed in Sec. 4.2 and shown 



in Fig. 11. 



The energy and momenta of the ?;-sector particles not decaying back into the SM is 
the prime source of the missing p±_ distributions, ^ see Fig. ^. In each of the six setups 
there is only one source of missing energy. In the Abelian ones it is the q^s that escape 



^with some extra effects from neutrinos e.g. in 6, c and r decays, included in the plots but here not 
considered on their own. 
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Figure 7: Mz' vs KM^^: the energy spectrum of the 7^ and diagonal ■n^l Pv emitted per event. On 
top: left side shows the energy distribution for AM^' , the right side shows the corresponding one 
for NAM^' . Bottom: left side shows the energy distribution for KMA-y^ , right side shows KMNA^^ . 
For the abelian cases mq^ = 50 GeV, for the non abelian cases rriq^ = m-^^/i. In all four cases 
Uv = 0.4. 



detection, while in the non- Abelian ones it is the stable non-diagonal v-mesons. For KMA^^ 
the falling Abelian spectra are easily understood from the bremsstrahlung nature of 
the 7^ emissions. The spike at = comes from events without any emissions at all, 
where all the energy is carried away by the invisible q^s, and would hardly be selected by a 
detector trigger. (ISR photons might be used as a trigger in this case, but with irreducible 
backgrounds e.g. from — )• uu it is not likely.). In the KMNA case the momentum of 
non-diagonal f-mesons does not leak back, this again allows a falling slope and a spike at 
/i_L = 0, for events in which equal amount of energy in the non-diagonal mesons radiated 
from either side of the q^q.^ system. For the SMA and SMNA scenarios, on the other hand, 
the starting point is the pT imbalance that comes from the and e~ from the and 

decays, which have no reason to balance each other. So even without 7^, emission, or 
diagonal tTv/ Pv^ there will will be a p± imbalance. 

In the SMNA case, though the spectrum is shifted towards lower missing px_ = 
because on average a higher number of mesons are radiated, so it is less likely to have an 
event with the two leptons back-to-back. There could also be = cases in which all 
the mesons are flavour diagonal and all the energy-momenta decays back into the SM, but 
these events are very rare. 
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Figure 8: KMAj^, KMNA^^, SMA, SMNA: the spectrum in each event. For the Abehan cases 
(left) the ■j^^ is due to the escaping detection. For the non-Abehan cases (right) it is due to the 
u-flavoured mesons not decaying into SM particles. In the Abelian cases ruq^ — 50 GeV, while in 
the non- Abelian cases (right) ruq^ = 1, 3, 5, 10. For all plots a„ = 0.4. 



In the Abelian case, the missing distribution is directly connected to the mass 
parameter values m^^ and, in the SMA case, to the niE^- The value of m^^ in the KM/ Z' 
mediated cases may be extracted from the kinematic limit given by the "shoulder" of 
the distribution. In the SM-mediated case, where two different fermion mass scales are 
involved, one can extract a relationship for the relative size of the two from lepton energy 



distributions such as the one in Fig. 16, see for details. The distribution that directly 
pinpoints the mass of the particle decaying back into the SM, though, is the invariant 
mass of the lepton pairs produced, and that of the hadronic jets. We will discuss these 
distributions in the sections dedicated to each scenario. 

4.2 AM.Z' and NAM^' 

In discussing the phenomenology of the different scenarios we will describe the f-sector 
particle distributions first, then the visible particle distributions followed by the jet distri- 
butions. 

The number of particles of 7^, photons emitted in the AM^/ and NAM^/ scenarios 
was described in Fig. Q and Fig. ^ in the previous section. The difference between the 
Abelian and non-Abelian dependence on the and fn^^j^^jp^ parameters has already 



been highlighted in the same Sec. |4^, as well as the ^v/t^v/ Pv energy distributions, the 
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Figure 9: AMz' and NAMz'- The distribution of the invariant mass of the lepton pairs. Note the 
peak at 2 GeV, in both cases corresponding to the mass to be reconstructed to^„ — nip^/^^. On 
the left, in the Abelian case, niq^, = 50 GeV, while on the right, in the non-Abelian case, m^^ = 1 
GeV. The coupling is fixed at a„ = 0.4. 



charged multiplicity and the spectrum. The difference between KMA^^ and SMA was 
also discussed. 

The difference between the Abelian and non-Abelian distribution in Fig. ^ is more 
subtle. In the Abelian case an event has maximum p± inbalance when one of the Qv/Qv 
produced emits a collinear 7^ which takes most of the momentum while the other 

v-quark has no emission and goes undetected. The more 7^, are emitted, the less likely 
it is that the undetected will have maximal energy. This remains true for all the m^^ 
contemplated (except in the low-p_L region). In the non-Abelian case, to have large p± 
inbalance the event must produce few energetic mesons back-to-back and have the mostly 
flavoured mesons at one end and mostly flavour neutral mesons at the other end. When the 
meson mass is lower, there is a higher probability of the string producing a large number 
of mesons and the likelihood of having large fi± falls rapidly. When the meson masses are 
higher and fewer TT^/p^ are produced the high ^_|_ distribution falls off less rapidly. 

The 'jv/'^v/Pv mass can be extracted from the lepton pair invariant mass, where it 
shows up as a well-defined spike, Fig. (The additional spike near zero mass is mainly 
related to Dalitz decay vr'' — )• e~^e~'y.) Once the mass is known, the remaining hadrons and 
photons may be clustered using the Jade algorithm, with m^^/j^^/p^ as the joining scale. 
The corresponding number of jets and invariant mass distribution for the hadronic jets is 
given in Fig. The jet invariant mass distribution clearly shows the peaks connected 
to the invariant mass of the 7^. The background comes from several sources. The spike 
at zero mass is mainly related to ISR photons; although we assume no detection within 
50 mrad of the beam directions, some isolated photons do show up above this angle and 
form jets on their own. When kinematically possible, r decays and c and b decays will also 
occur. These contain neutrino products that reduce the visible mass, thus contributing to 
a continuum below the mass peak. Finally, misidentifications among partly overlapping 
systems leads to tails on both sides of the peak. 

An efficient clustering algorithm should maintain the ratio between the average number 
of 7t, particles produced and the number of jets found, as is confirmed by comparing the 
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Figure 10: AMz'- the number of jets per event and the distribution of the invariant mass of the 
of the jets. In both plots the mass is m^^ = 50 GeV. The left side shows the dependence 
for rrij^ = 10 GeV, while the right side shows the jet invariant mass distribution for = 0.4 and 
TO^^ — 6 GeV. The jet reconstruction algorithm is Jade, with rricnt corresponding to the 7„ mass. 
Note the peak at 6 GeV. 



plots on left side of Figs. ^ and [l^. In this particular exercise we have relied on the 
extraction of the relevant mass scale from the lepton pair invariant mass distribution. This 
should be guaranteed by the presence of leptons in all six scenarios. Specifically, if the only 
way to decay back into the SM is via kinetic mixing, the 7^ — )• SM branching ratios are 
fixed by the off-shell 7* branching ratios. In the non-Abelian case the absence of a spin 1 
i;-meson would reduce the rate of e~^e~ and /U^^^ pairs by helicity suppression. However, 
also in such scenarios, a simple trial-and-error approach with a range of jet clustering scales 
would suffice to reveal a convincing jet mass peak. 

More information about the event and the model parameters can be extracted from 



the angular distributions in Sec. 4.5. 



4.3 KMA^„ and KMNA^^ 

The Z' mediated and 7^ mediated setups are effectively very similar, once the difference 
in coupling constants is factored out. The phenomenology can appear rather different, 
however, because of the initial state radiation from the electron/positron beams. To view 
this, recall that the photon bremsstrahlung spectrum is spiked at small energy fractions, 
oc dz^/z^, and that therefore the electron- inside-electron PDF roughly goes like dze/ (1— Ze), 
with Ze = 1 — z^. In the case of a 7 or light propagator, behaving like 1/s = l/{zes), this 
combines to give a dze/ (^e(l — ^e)) = dz^/ ((1 — z^)z.y) spectrum. The complete description 
includes the emission of multiple photons off both incoming beams, but the key features 



above described are readily visible in Fig. 11. Specifically, the spike at -Eisr = 400 GeV 



corresponds to the emission of an energetic photon on one side only, while the non-negligible 
tail above that requires hard emissions on both sides. 

Note that while the number of u-particles produced are similar for KMA.y^ and AM^/ 
in the Abelian case. Fig. ^ vs. Fig. |l^, they are different in the non-Abelian case. Fig. ^ 



vs. Fig. 13. Specifically, the higher average multiplicity in the non-Abelian case leads to 
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Figure 11: KMNA^^: The total energy radiated in ISR photons. Note the spike around 400 GeV: 
a large fraction of the events will have a reduced s due IS photon emission. The mass is fixed 
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Figure 12: KMA^^: the number oi^^ gauge bosons emitted per event. On the left we emphasize 
the ckt, dependence, while on the right the m-y^ dependence. In both plots the mass is niq^ = 10 
GeV. On the left side m^^ = 10 GeV and = 0.1, 0.4, 0.6, while on the right side m^^ = 2, 6, 10, 20 
GeV and the coupling is fixed at = 0.4. 



a double spike in the distribution, a normal one from events with little ISR and an extra 
low-multiplicity one from events with much ISR. 

4.4 SMA and SMNA 

The number of f-particles produced in the standard model mediated scenarios is shown 
in Fig. 14, and the energy of the 7t, photons and of the u-mesons in Fig. The kinetic 



boundary E^^^ is different in the two cases, owing to the choice of niq^ = 50 GeV in the 
Abelian case. This reduces the energy available for 7„ emissions. The |f_L distribution was 
shown in Fig. |^ and has already been discussed. 

The most important distribution to pinpoint the masses of the model is the lepton 
energy spectrum, Fig. In this case leptons may come from both the kinetic mixing 
decays of the 7^ and from the decays of the into eq^j. The energy spectra are very 
different in the two cases. The leptons coming from the E^ decays tend be highly energetic. 
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Figure 13: KMNA^^: the number of flavour diagonal -KvI Pv gauge bosons emitted per event. On 
the left we emphasize the dependence, while on the right the Tn-j^^/p^ dependence. On the left 
side m.;r„/p„ = 10 GeV and Uy = 0.1, 0.4, 0.6, while on the right side m^^/p^ = 2, 6, 10, 20 GeV and 
the coupling is fixed at — 0.4. 
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Figure 14: SMA: the number of gauge bosons emitted per event. On the left we emphasize 
the a^j dependence, while on the right the m^^ dependence. In both plots ttie^ = 250 GeV and 
niq^ = 50 GeV. On the left side m-y^ = 10 GeV and Uv = 0.1,0.4,0.6, while on the right side 
m-y„ = 2, 6, 10, 20 GeV and the coupling is fixed at a„ = 0.4. 



while the rest are less so. A reasonable first approximation is to associate the highest energy 
electron and positron with the two decays, and the rest with the Jv/t^v/Pv ones. The 



curve in Fig. 1(: represents the sum of the steeply falling spectrum associated to the leptons 
coming from 7„ decay, and a flat spectrum associated to the leptons from — )• eq^ decay. 
The upper and lower shoulders of the former energy distributions then give a relationship 



between the E^ and the masses |18]. 



4.5 Angular distributions and event shapes 

The distribution of the production cross-section as a function of the angle between the 
jets and the beam axis has a characteristic dependence on the spin of the pair-produced 
particles. This fact may be used to identify the spin. 



In Fig. 17 one may observe the cos 9i distribution, where 9i is the polar angle between 
the i^^ jet and the beam axis, for m^^ = 10 GeV. Only jets with a reconstructed mass 
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Figure 15: SMA vs SMNA: the energy spectra of the 7^, and of the diagonal ttv/ Pv emitted per 
event. The left side shows the energy distribution for SMA m^^ = 250, GeV m^^ = 50 GeV, 
the right side shows the corresponding distribution for SMNA (in this case ruq^ = 1, 3, 5, 10 GeV 
corresponds to half of the i^vl Pv mass.). 
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Figure 16: SMA vs SMNA: the lepton energy spectra in the Abelian (left) case and in the non- 
Abelian (right) case. In the Abelian case GeV — 50 GeV, while in the non-Abelian case 
ruq^ = 1, 3, 5, 10 GeV. In both cases ruE^ = 250 GeV. 

rrij > are shown, since lower-mass jets are strongly contaminated by ISR photons 

above the 0cut = 50 mrad cut. The production cross-section for e~^e~ — )• qvqv, with g„ a 
massless spin 1/2 fermion, is proportional to 1+cos^ 9. In the AMz' case one must allow for 
corrections due to the 6'cut; to the being massive and to the 7„ radiation; typically this 
leads to a somewhat flatter distribution. In the SMA case the isotropic decays — t- eq^ 
flattens whatever original e'^e~ — )■ distribution, unless the eqy decay products are 

highly boosted. In our case a small ISR contamination is still visible close to cos^ = 1, 
but otherwise the distribution is flat. 

The distribution of opening angles in pairs of u-particles, 9ij, should give some insight 
whether the secluded sector G is an Abelian or a confining non-Abelian group. In an 
Abelian event the q^q^ quarks define a dipole emission axis. To first approximation the 
t;-gammas are emitted independently, i.e. with a flat distribution in the (p angle around 
the dipole axis, and uniformly in rapidity along this axis. In the non-Abelian case the 
emissions occur along a chain of dipoles, that is reconfigured by each new emission, since 
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Figure 17: AMz' and SMA: the distribution of cos6'i, where 9i is the angle between the i^^ jet 
and the beam axis, for m^^ = 30 GeV and for niq^ — 50 GeV. A lower cut on the invariant mass 
rUj > m.^^/jr„/p„ /2 was applied in this case. In both cases m.^^, = 10 GeV and the u-coupling is 
fixed to a„ — 0.4, iriq^ — 10 GeV. 
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Figure 18: AAf^ and NAAf^: the distribution of 6ij, the angle between the 7^s in the Abelian case 
(left), or the angle between the diagonal mesons in the non- Abelian case (right). The distributions 
are shown as functions of the mass m-y^ = itIt^^/p^ . The coupling constant is = 0.4 in both cases. 
Note that the number of and diagonal it„/ is different in the two cases. 



the g^s carry u-colour charge. This implies a different underlying correlation structure, but 
it is unclear what happens with this correlation on the way through the v-hadronization 
process and the decays back into the standard sector, and how best to search for it. 

In Fig. 18 we show the Oij distribution of u-particle pairs, for the Abelian AM^/ and 
non- Abelian NAM^/ cases. Note how the 9^ and Oij ir angles are preferred in non- Abelian 
case. The comparison is somewhat misleading, however, in the sense that we compare 
scenarios with the same niy^ = m^^/p^ and the same ay, but with different numbers of 
t;-particles per event (see Fig. ^ and different u-particle energy distributions. In Sec. ^ we 
will discuss further these distributions under more similar conditions. 

In Fig. |l^ we show the reconstructed jet-jet cos 6ij distributions corresponding to the 
?;-particle distributions in Fig. 18. Note how the relative difference between the Abelian 
and non- Abelian scenarios is maintained. 

In order to characterize the shape of the events one may also use thrust and (the 
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Figure 19: AM^ and NAM'^: the distribution of 9ij, the angle between the jets in the Abelian 
case (left) and in the non- Abelian case (right). The distributions are shown as functions of the 
mass rrij^ = m^^/p^. The u-coupling is fixed to = 0.4. 



linearized version of) sphericity [39, 4C, 41 1. These indicate whether an event is more pencil- 
like, T=l and S = 0, or more spherical, T = 1/2 and S = 1. Sphericity and thrust are 
primarily intended for events analyzed in their own rest frame, while the visible systems we 
study have a net momentum that is compensated by the stable secluded-section particles, 
plus ISR photons going down the beam pipe and neutrinos. Since we are interested in the 
properties of the visible system itself, not in its net motion, the analysis is performed in the 
rest frame of this visible system. In addition, for the SMA scenarios, a further distortion 
occurs by the kinematics of the — )• eq^ decays, and by the presence of the resulting 
in the final state. To this end, the highest-energy electron and positron are excluded from 
the definition of the visible system. 

In Fig. we show the thrust distributions in the six scenarios, and in Fig. 21 the 
sphericity distributions. As one may have predicted, the events become less pencil-like as 
the coupling grows. In addition the SM mediated events are less likely to be pencil-like 
than the KM or Z' mediated ones. Note that events with smaller Oy, in which nothing is 
radiated are not analyzed. 



5. Analysis: comparing ^(1) and SU{N) 

In this section we begin to address the issue of discriminating between Abelian and non- 
Abelian scenarios in cases in which smoking- gun discriminating signals are absent. To this 
end we consider the most challenging scenario, in which and ir^, have the same mass, 
and the same average number of u-particles leak back into the SM sector, carrying the same 
average amount of energy. We also consider the same production mechanism, to reduce 
model dependence and to isolate the effects of the hidden sector dynamics. 

So how could the Abelian and non-Abelian scenarios come to have so similar proper- 
ties? First off, of course, = m^r^ = rup^ needs to be assumed. Thereafter the value of 
A'fiav in the non-Abelian model specifies that exactly an average fraction l/A'gav of the full 
energy leaks back into the visible sector. In the Abelian model, for a given m^^, Oy is the 
only free parameter that could be fixed to give that average energy. For this parameter 
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Figure 20: The plots show the thrust distribution for the six scenarios as a function of the a^- 
These correspond to, from top to bottom, Mz', KM7„ and SM production. The plots on the left 
show the Abelian case while the plot on the right show the non-Abelian cases. For all Abelian plots, 
the parameters are set to uig^ = 50 GeV and m^^ — 10 GeV. For the SMA case ttie^, = 250 GeV is 
set as well. For the non-Abelian cases, the mass of the mesons is fixed to fn^^^/j^^/p^ 



set, the number and energy spectrum of 7^s are predicted entirely by the perturbative cas- 
cade. These distributions now need to be roughly reproduced by the non-Abelian model, 
which first of all means the same average number of f-particles decaying back into the 
SM. While nig^ = is fixed in this case, there is freedom in the choices of and 

non-perturbative fragmentation parameters. Recall that the number of v-particles will not 
vanish in the ^ limit for the non-Abelian model, unlike the Abelian one. Actually it 
turns out to be slightly difficult to reduce the non-Abelian multiplicity down to the level 
set by the Abelian scenario. With an comparable to that of a QCD cascade at a cor- 
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Figure 21: The plots show the (hnearized) sphericity distribution for the six scenarios as a function 
of the ay. These correspond to, from top to bottom, KM^/, KMjy and SM production. The plots 
on the left show the Abelian case while the plot on the right show the non-Abelian cases. For all 
Abelian plots, the parameters are set to m^^ = 50 GeV and m-y^, = 10 GeV. For the SMA case 
niE^ = 250 GeV is set as well. For the non-Abelian cases, the mass of the mesons is fixed to 



responding energy/mass ratio, the longitudinal fragmentation function needs to be made 
harder by decreasing a and increasing h' relative to the QCD values. 

Using such a strategy, a matching pair of scenarios have been constructed, an AM^ 
model with ruq^ = 20 GeV, m^^ = 10 GeV and a„ = 0.3, and a NAM'^ with niq^ = 5 GeV, 
^TTv/pv ~ GeV, A^flav = 4, Ot, = 0.15, a = 0.12 and h' = 2. This gives fair agreement, as 
can be seen in Fig. p2| . 

We can now compare the angular distributions in the two cases. In Fig. ^ (left) one 



z 

* 



0.4 
0.35 

0.3 
0.25 

0.2 
0.15 

0.1 
0.05 




Ab 
non-Ab. 



10 



1 

0.1 
0.01 
0.001 
0.0001 
le-05 









Ab. ' — — <' 




non-Ab. " ' 
















I 



100 200 300 400 
E, [GeV] 



500 



u 
Z 
* 





100 200 300 400 500 
pT [GeV] 



Figure 22: The plots show the eomparison between the Abehan and the non-Abelian setups: 
the number of j^s (Abehan) or diagonal iTys/pyS (non-Abelian) (top left), the w-particle energy 
distribution (top right), the number of SM charged particles produced, and the missing p±. The 
scenarios are chosen to yield similar distributions. 
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Figure 23: The plots show the comparison between the Abelian and the non Abelian setups. On 
the left is the cos 9ij between the w-particles which can decay back into the SM. On the right is the 
corresponding cos 9ij between the jets in the detector. 



may observe how in the non-Abelian case the 9ij distribution of the angle between the 
visible f-particles is much more peaked near and vr than in the Abelian one. The plot on 
the right shows how this characteristic is maintained in jet distributions. 

We have repeated the study for some different model parameters and found sumilar 
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results. In the seeond comparison the AM^ model has niq^ = 50 GeV, = 10 GeV 
and ay = 0.43, while NAM'^ has ruq^ = 5 GeV, rrij^^/rhov — 10 GeV, A^flav = 4, = 0.18, 
a = 0.2 and 6' = 2. 

It is tempting to ascribe the observed differences to different radiation/hadronization 
patterns in the two scenarios, ultimately deriving from the different dipole emission topolo- 
gies already discussed. One attractive possibility is that the relative lack of peaking at 
cos 9ij = ±1 for the Abelian scenario is a consequence of the dead cone effect, i.e. the 
suppression of emissions parallel to a massive radiating particle. A back-to-back q^q^ 
pair undergoing non-Abelian hadronization would have no corresponding suppression for 
?;-hadron formation along this axis. However, considering of the flexibility that exists in 
the tuning of non-perturbative hadronization parameters, and the differences observed in 
the i;-particle energy spectrum, we will not now be as bold as to exclude the possibility 
of a closer match. If such a match required straining the non-perturbative model to be- 
have rather differently from QCD extrapolations, however, then such tunes would not be 
particularly credible. 

6. Summary and Outlook 

In this article we have compared six different scenarios of a generic secluded-sector char- 
acter. These are either kinetic mixing with a light 7^, or Z' , or new Fy particles, and in 
each case either with a broken U{1) or an unbroken SU{N). The Fy particles are charged 
both under the Standard Model groups and the new secluded-sector groups, and so are 
guaranteed a significant production rate, whenever kinematically possible, whereas the rate 
via '-fy or Z' depends on a number of parameters such as the 7/71, mixing parameters or 
the Z' coupling structure and mass. In this article we have completely disregarded such 
rate issues and instead studied the properties of the different models on an event-by-event 
basis. 

In order to do so, we have developed a new flexible framework that implements 
hadronization in the hidden sector. Similar modeling in the past have relied on simple 
rescaling of QCD, whereas here we set up hidden-sector string fragmentation as a com- 
pletely separate framework, though sharing the same underlying space-time structure of 
the hadronization process. We have also expanded on our previous work with parton show- 
ers in the hidden sector, possibly interleaved with radiation in the visible sector, by allowing 
for the emission of massive 7^s, when the U{1) group is broken. In order to obtain the 
correct behaviour, both in the soft/collinear limits and for hard emissions, the shower is 
matched to first-order matrix elements we have calculated for the massive 7^, cases. Much 
of the framework presented here could be applied also to other related scenarios, although 
there are limits. For instance, implicitly it has been assumed that the Qy masses are not too 
dissimilar from the confinement scale A of the new SU (N) group — hadronization would 
look rather different in the limit niq^ S> A. 

In the scenarios where FyS are produced and promptly decay like Fy — )■ fqy, the 
presence of the fermions / in the final state is a distinguishing factor, and the / energy 
spectrum can be used to extract information on mass scales in the secluded sector. At 
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first glance, the production mediated by a 7^, or a Z' would seem to be more similar. 
The different location of the propagator mass peaks leads to quite significant patterns of 
initial-state photon radiation, however, that would be easily observed. 

The challenge, thus, is to distinguish an Abelian and a non-Abelian scenario interac- 
tions in the secluded sector. In certain cases that would be straightforward, e.g. if there 
is only one g^, species, so that all energy decays back into the standard sector in the non- 
Abelian models. It is possible to fix parameters in the non-Abelian cases so that only some 
fraction decays, so that the level of activity matches Abelian one. 

Our first investigations here point to differences showing up in event properties related 
to the overall structure of the energy and particle flow. More elaborate tunings possibly 
might bring the models closer together, but one would hope that data still would favour 
"what comes naturally" in cither of the models. 

There are also more handles than the ones we have used. We have shown that not only 
lepton pairs but also jets are amenable to mass peak identification, which would allow to 
divide events into several subsystems and analyze their relative location, e.g. searching for 
coherence effects. Should lifetimes be long enough to induce displaced vertices, not only 
would that be a spectacular signal, but it would also be a boon to such analysis efforts. 

The most obvious next step would be to study these models for consequences at the 
LHC. The task can be split into three parts: cross sections, triggers and model-specific 
event properties. The cross sections are so intimately related to the choice of masses and 
couplings that it will be impossible to exclude the possibility of a secluded sector, only 
to separate excluded and not-excluded regions of parameter space, in close analogy with 
SUSY. The obvious trigger would be but we have seen that this would not work 
for non-Abelian scenarios with one g^, flavour. It would then need to be supplemented 
by the presence of (multiple) lepton pairs of some fixed invariant mass and, if we are 
lucky, displaced vertices. The final step would be to understand whether the more busy 
environment in hadronic events would still allow to distinguish Abelian and non-Abelian 
models — the separation between the three production scenarios we have considered here 
should still be straightforward. Chances are that this will bring us full circle to the cross 
sections issue, since more sophisticated analyses will require a decent event rate to start 
out from. 
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A. Scencirio selection and setup 

We here present some information relevant to get going with secluded-sector event gener- 
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name 


name 


identity 


comment 




Dv 


4900001 


partner to the a quark 


(A) 


Uv 


yi r\ r\ r\ r\ r\ r\ 

4900002 


partner to the u quark 


Sv 


Sv 


4900003 


partner to the s quark 




Cv 


4900004 


partner to the c quark 




Bv 


4900005 


j_ j_ j_1 7 1 

partner to the o quark 


Tv 


Tv 


4900006 


partner to the t quark 




Ev 


4900011 


partner to the e lepton 




nuEv 


4900012 


partner to the neutrino 


My 


MUv 


4900013 


partner to the iepton 




nuMUv 


4900014 


partner to the neutrino 


Tv 


TAUv 


4900015 


partner to the r iepton 




nuTAUv 


A r\r\r\r\ ■'K 

4900016 


partner to the neutrino 


9v 


gv 


4900021 


the f;-giuon m an dc/ (iV j scenario 


Iv 


gammav 


4900022 


the v-plioton in a scenario 




Zv 




massive gauge boson linking SM- and w-sectors 


Iv 


qv 




matter particles purely in f-sector 


diag 


pivDiag 


4900111 


flavour-diagonal spin i;-meson 




rhovDiag 


4900113 


flavour-diagonal spin 1 t!-meson 


up 


pivUp 


4900211 


flavour-nondiagonal spin f-meson 


up 
Pv 


rhovUp 


4900213 


flavour-nondiagonal spin 1 f-meson 






4900991 


glueball made of u-gluons 



Table 2: The allowed new particles in valley scenarios. Names are gives as used in this text and 
as in Pythia 8 event listings. The identity code is an integer identifier, in the spirit of the PDG 
codes, but is not part of the current PDG standard 

ation in Pythia 8. Basic knowledge of the program is assumed |25| ]. 

The f -particle content is summarized in Tab. |2[ Needless to say, not all of them would 
be relevant for each specific scenario. Internally further copies of may be used, up to 
code 4900108, but these do not appear in the event record. Properties of the particles can 
be set to modify the scenarios, notably mass (mO); only the must remain massless. If 
Fv — ^ f<lv is allowed, masses must be chosen so that the decay is kinematically possible. 
The TT^ and p^ masses should be set at around twice the one. For unstable particles the 
width (mWidth) and allowed mass range (mMin and mMax) can be set. To generate displaced 
vertices the cr value must be set (tauO). Spin choices are described later. 

Several particles by default are set stable, so it is necessary to switch on their decay 
(mayDecay). For 7^, vrj^"^^ and p^'^'^ no decay channels are on by default, since that set 
depends so strongly on the mass scale selected. The simple way, of switching on everything 
(onMode = on) works in principle, but e.g. a 10 GeV 7^, would then be above the hb 
threshold but below the BB one, and so generate a trail of (harmless) error messages. Also 
the branching ratios of the decay channels may need to be adjusted, based on the scenario. 
The Zy ones are mainly place-fillers, to give one example. 
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1 

code 


flag 


process 


4901 


HiddenValley 


gg2DvDvbar 


99 - 




A nno 


HiddenValley 


gg2UvUvbar 


99 - 


T T TT 




A nno 

491)3 


HiddenValley 


gg2SvSvbar 


99 - 


V c c 




49(J4 


HiddenValley 


gg2CvCvbar 


99 - 






49Uo 


HiddenValley 


gg2BvBvbar 


99 - 






49Ud 


HiddenValley 


gg2TvTvbar 


99 - 






49ii 


HiddenValley 


qqbar2DvDvbar 


qq- 






4912 


HiddenValley 


qqb ar 2UvUvb ar 


qq- 


> 9* 




4913 


HiddenValley 


qqbar2SvSvbar 


qq- 


>5* ^ 


SvSy 


A r\i A 

4914 


HiddenValley 


qqbar2CvCvbar 


qq- 


>5* ^ 


CvCy 


4915 


HiddenValley 


qqbar2BvBvbar 


qq- 


>5* ^ 


By By 


4916 


HiddenValley 


qqb ar 2 T vT vb ar 


qq- 


>5* ^ 




4921 


HiddenValley 


ffbar2DvDvbar 


ff_- 


-> 7* — ) 


DyDy 


4922 


HiddenValley 


f fbar2UvUvbar 


ff- 


->7* ^ 


UyUy 


4923 


HiddenValley 


ffbar2SvSvbar 


ff- 


7* — ) 


SySy 


4924 


HiddenValley 


ffbar2CvCvbar 


ff- 


7* — ) 


CyCy 


4925 


HiddenValley 


ffbar2BvBvbar 


ff- 


->7* ^ 


By By 


4926 


HiddenValley 


ffbar2TvTvbar 


ff- 


* 

->7 ^ 




4931 


HiddenValley 


ffbar2EvEvbar 


ff - 


* 

->7 ^ 


EyEy 


4932 


HiddenValley 


f f bar2nuEvnuEvbar 


ff- 


* 

->7 ^ 




4933 


HiddenValley 


ffbar2MUvMUvbar 


ff- 


->7*^ 


My My 


4934 


HiddenValley 


f f b ar 2nuMUvnuMU vb ar 


ff- 


->7*^ 


I' My i^My 


4935 


HiddenValley 


ffbar2TAUvTAUvbar 


ff- 


7* — ) 


TyTy 


4936 


HiddenValley 


ffbar2nuTAUvnuTAUvbar 


ff- 


7* — ) 




4941 


HiddenValley 


ffbar2Zv 


ff- 


^:(- 





Table 3: Allowed processes that can be switched on individually. 



The list of processes is shown in Tab. |^. It would be possible to switch on all of them 
with HiddenValley : all = on, but normally that would imply a mix of different scenarios 
that does not appear plausible. Many processes should also be viewed in the context of 
the other choices made. 

Finally, the list of relevant model parameters is shown in Tab. ^ On top is the 
choice between a U{1) and an SU{N) scenario. The Fy and qy spins must be selected 
in a coordinated fashion, to be consistent with Fy — )• fqy decays. The choice of Fy spin 
directly affects the process (differential) cross sections. If Fy has spin 1 also the choice of 
an anomalous magnetic moment would have an influence. 

The kinetic mixing switch allows to reuse the 7*-mediated processes in a completely 
different context than originally foreseen, in which the Fy have no Standard Model coupling 
but are produced by 7 — 7^, mixing. Actually this redefines the Fy to be equivalent with 
what we normally call qy. Thus a process like ff^J*^ EyEy becomes ff^^*^ 
it ~^ QvQv To complete this transformation you need to set the Ey stable (mayDecay = 
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parameter 


def. 


meaning 


Scenario 


HiddenValley : Ngauge 
HiddenValley : spinFv 
HiddenValley : spinqv 
HiddenValley : kappa 
HiddenValley : doKinMix 
HiddenValley : kinMix 


3 

1 



1. 

off 

1. 


1 for C/(l), N for SU{N) 

0, 1 or 2 for Fy spin 0, 1/2 and 1 
spin or 1 when sp^ = 1/2 
anomalous magnetic dipole moment 

allow kinetic mixing 

strength of kinetic mixing, if on 


Showers in secluded sector 


HiddenValley :FSR 
HiddenValley : alphaFSR 
HiddenValley : pTminFSR 


off 
0.1 
0.4 


allow final-state radiation 
constant coupling strength 
lower cutoff of shower evolution 


Hadronization in secluded sector 


HiddenValley : fragment 
HiddenValley : nFlav 
HiddenValley : probVector 
HiddenValley : aLund 
HiddenValley : bmqv2 
HiddenValley : rFactqv 
HiddenValley : sigmamqv 


off 
1 

0.75 

0.3 

0.8 

1.0 

0.5 


allow hadronization 

A^fiavi number of distinct q^j species 

fraction of spin-1 f -mesons 

a parameter in eq. (A.l) 

b' = bni'^^ parameter in eq. (A.l) 

r parameter in eq. (A.l) 

a', such that a = cr'niq^ 



Table 4: The parameters that can be set to select the model to study, with default values and 
some expanations. 



false), uncharged (chargeType = O) and invisible (isVisible = false). 

The shower parameters should be self-explanatory. The lower cutoff scale can be raised 
in proportion to the characteristic mass scales, but ought to be no more than mq^/2, say. A 
lower cutoff means longer execution time without any significant change of event properties. 

The hadronization parameters have also been discussed before, except for r, which 
is providing slightly more flexibility to the Lund-Bowler fragmetation function than in 
eq. ( p6D 

/(^) « ^ (1 - exp ^ . (A.l) 

where r could be set anywhere between and 1. The dimensionless a' parameter is normal- 
ized so that the q^, of each new pair produced in the hadronization has a {p']_) = (cr'mg^)^. 

Behind the scenes, the HiddenValleyFragmentat ion class can reuse most of the stan- 
dard StringFragmentation and MiniStringFragmentation machineries. Specifically, al- 
ready for the Standard Model hadronization, the selection of flavour, z and p± is relegated 
to three "helper" classes. The three new classes HVStringFlav, HVStringZ and HVStringPT 
derive from their respective SM equivalent, and cleanly replace these three aspects while 
keeping the rest of the handling of complex string topologies. Specifically, it would be 
straightforward to expand towards a richer flavour structure in the secluded sector. Note, 
however, that it is important to select f-quark "constituent" masses that reasonably match 
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the intended ?;-meson mass spectrum, since such relations are assumed in parts of the code. 
Furthermore, with new defined with separate particle data, one must disable the few 
lines in HiddenValleyFragmentation: :init(. . .) that now duplicate the one into 
several identical copies. 
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